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We explore the phenomenon of spiral defect chaos in two types of generalized Swift-Hohenberg 
model equations that include the effects of long-range drift velocity or mean flow. We use spatially- 
extended domains and integrate the equations for very long times to study the pattern dynamics 
as the magnitude of the mean flow is varied. The magnitude of the mean flow is adjusted via a 
real and continuous parameter that accounts for the fluid boundary conditions on the horizontal 
surfaces in a convecting layer. For weak values of the mean flow we find that the patterns exhibit 
a slow coarsening to a state dominated by large and very slowly moving target defects. For strong 
enough mean flow we identify the existence of spatiotemporal chaos which is indicated by a positive 
leading order Lyapunov exponent. We compare the spatial features of the mean flow field with that 
of Rayleigh-Benard convection and quantify their differences in the neighborhood of spiral defects. 

PACS numbers: 47.54.-r, 47.52. +j, 05.45.Jn, 05.45.Pq 



I. INTRODUCTION 



The chaotic behavior of spatially-extended dissipative 
systems has been intensively studied in recent years [ij. 
Spatiotemporal chaos has been observed in a wide range 
of physical systems including Faraday waves that appear 
on the surface of an oscillating layer of fluid ^] , reacting 
chemical mixtures [3] , excitable media Q , and Rayleigh- 
Benard convection in a shallow fluid layer heated from 
below [ll]. In particular, the study of Rayleigh-Benard 
convection continues to provide fundamental insights into 
the dynamics of pattern forming systems that are driven 
far-from-equilibrium [l|, Q- The state of spiral defect 
chaos has received significant attention since its discovery 
by Morris et al. ^ . Spiral defect chaos is characterized by 
the complex dynamics of rotating spiral defects and inter- 
estingly, occurs for fluid parameters where straight par- 
allel convection rolls are linearly stable Jj, i5|] . Although 
much effort has been spent on building our understand- 
ing of the origins and dynamical features of spiral defect 
chaos, many open questions still remain. 

A signiflcant difficulty in studying spatiotemporal 
chaos is that the experimental systems of interest are 
often large and strongly driven. This presents significant 
obstacles for both analytical and numerical approaches. 
For example, the numerical simulation of the Boussinesq 
equations, which govern Rayleigh-Benard convection, in 
large domains, for long simulation times, and for many 
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values of the system parameters is out of reach using cur- 
rently available algorithms and computing resources. Al- 
though significant progress has been made in the ability 
to simulate convection for experimental conditions, the 
computational cost remains very high [Ti] . In light of chal- 
lenges such as these, the use of simpler model equations 
has played a pivotal role in furthering our physical un- 
derstanding of spatiotemporal chaos. In Rayleigh-Benard 
convection, the two-dimensional Swift-Hohenberg equa- 
tion [8| has led to numerous physical insights regarding 
questions of pattern formation [l|. However, the use of 
the Swift-Hohenberg equation to study spatiotemporal 
chaos, and in particular spiral defect chaos 3, has been 
called into question [10], as will be further discussed be- 
low. This leaves no clear choice for a model system to 
be used for the study of spiral defect chaos with direct 
relevance to fluid convection. 

The presence of a long-range mean flow is well known 
to play an important role in the dynamics of Rayleigh- 
Benard convection ^, JJj, JJi] . This led to extensions of 
the Swift-Hohenberg equation to account for these ef- 
fects [3, |lj, ^^. The mean flow is a weak horizontal 
flow field that acts over length scales larger than that of 
the convection rolls; it results from the coupling to fluid 
vertical vorticity and is induced by roll curvature, am- 
plitude gradients, and wavenumber gradients |ll. KLSl [l6| . 
The magnitude of the mean flow is much smaller than 
that of the convective roll motion, making it very diffi- 
cult to measure in experiment. A typical way to visualize 
the mean fiow in numerical work is to present contours 
of the integrated horizontal components of the fluid ve- 
locity over the depth of the fluid layer Q . Although the 
mean flow is weak it can have significant effect upon the 



dynamics and stability boundaries of the flow field and 
also adds a slow time scale to the dynamics. It has been 
shown numerically that the mean flow is required for spi- 
ral defect chaos in Rayleigh-Benard convection 17]. 

Spiral defect chaos was initially explored using numeri- 
cal simulations of the generalized Swift-Hohenberg equa- 
tion by Xi et al. [9|. These simulations were for rather 
short intervals of time t ~ 900 where t is the nondi- 
mensional time. We note that to generate the spiral de- 
fect state in this model system, the system needs to be 
taken far enough from the convective threshold; other- 
wise roll-type patterns would dominate with the system 
dynamics governed by the motion of topological defects 
such as grain boundaries, dislocations, and disclinations 
[18]. Schmitz et al. 10] later explored the generalized 
Swift-Hohenberg equation for the same parameters but 
for much longer simulation times t ^ 64000. Their re- 
sults suggested that the spiral defect state was only a 
transient with the long-time dynamics characterized by 
a slow coarsening process to a state dominated by large 
spirals. This work has cast doubt upon the ability of the 
Swift-Hohenberg equation to exhibit persistent dynamics 
that resemble spiral defect chaos. 

In this paper, we present a careful numerical study of 
these questions. In particular, we have performed very 
long-time simulations (t — 10^ time units) in very large 
domains, and for a wide range of mean flow strengths. 
We also compute the leading order Lyapunov exponent 
to determine if the dynamics are chaotic. Our investiga- 
tion is driven by the results on Rayleigh-Benard convec- 
tion that demonstrate the importance of mean flow upon 
spiral defect chaos 0, 0] . 

The remainder of this paper is organized as follows. 
In Sec. [H] we discuss two types of generalized Swift- 
Hohenberg model equations and provide some details re- 
garding its numerical solution and the computation of the 
leading order Lyapunov exponent. In Sec. lIIII we present 
our results and discuss the role of the mean flow upon 
the dynamics. Concluding remarks are given in Sec. IIVI 



II. APPROACH 
A. The Generalized Swrift-Hohenberg Models 

A generalized, dimensionless form of the two- 
dimensional time-dependent Swift-Hohenberg model is 
given by Q 

^+U-VV' = ei/'-(V2 + l)2v^ + Af[V'], (1) 



at 



a (V^ ~ c2 



V^C = g,n [V (VV) X Vij] ■ z, (2) 



where ip{x, y, t) is a scalar field describing the spatial and 
temporal variation of the convection patterns, e is a con- 
trol parameter giving the dimensionless distance from the 
convective threshold, and z represents the unit vector in 



the out-of-plane direction. The variable C(i, j/, t) is the 
vertical vorticity potential defined via $7^ = (V x U)^ = 
— V^Ci where Q,z is the vertical component of fluid vor- 
ticity and U is the mean flow or drift velocity. C, can also 
be interpreted as the stream function for the mean flow, 
given that 



u = V X (Cz) = (dyC, -d^o . 



(3) 



In Eq. ([2]) cr is proportional to the Prandtl number, g^ 
is a positive real constant that characterizes the mean 
flow coupling strength, and c is a real parameter intro- 
duced for modeling the effect of free-slip (c = 0) or no-slip 
(c ^ 0) boundary conditions for the horizontal surfaces 
of a convection layer. The term A^[^] in Eq. ([T]) rep- 
resents the nonlinearity of the system, which has many 
different forms in the literature based on system con- 
ditions [l|, [ij, [i9[. In this paper we present numerical 
results for two representative cases: (1) A^['0] = —"0^ 
for which the model equations are referred to as the 
Generalized Swift-Hohenberg (GSH) model [1^, and (2) 
N[il}\ = — (|V'0pV' + "0^) for which the equations are 
referred to as Manneville's model that was derived in 
Refs. [IJ, [lal from the Boussinesq equations. 



B. Numerical Simulations 

We numerically integrate Eqs. ffl-© using the ap- 
proach discussed by Cross et al. [2^ and provide only 
the essential details here. The domain is a square geom- 
etry that is discretized on a spatially uniform grid with 
periodic lateral boundary conditions. In our numerical 
simulations we begin from random initial conditions for 
-0 and set initially C,{x,y,t = 0) = 0. We discretize the 
spatial domain using a 512 x 512 grid with a grid spac- 
ing of Aa; = Ay = Ao/8 where Aq = 27r is the critical 
wavelength of a convection roll. These parameters are 
approximately equivalent to a Rayleigh-Benard system 
in a box geometry with an aspect ratio of F = L/d = 128 
where L is the side length of the box and d is the depth 
of the convection layer. Each individual simulation is al- 
lowed to evolve for t = 10^ time units using a time step 
of Ai = 0.2. 

The equations are stiff due to the very fast dynamics 
of the biharmonic operator in comparison to the much 
slower convective time scale of the pattern dynamics. 
An efficient solution is obtained using a pseudospectral 
operator-splitting approach |21| . The linear terms are 
treated exactly using an explicit exponential time in- 
tegration [221 l23j and the nonlinear terms are evolved 
forward in time using an explicit predictor-corrector ap- 
proach. To reduce the contributions of high wavenumber 
modes in the vorticity field a Gaussian filtering operator 
F-y is applied to the right-hand side of Eq. ([2]) [1^. In 
Fourier space it is given by F^ = exp f— 7^q^/2) where 7 
is the filtering radius and q is the wavenumber. In our 
simulations we have used a filtering radius of 7 = Ao/2. 



We compute the leading order Lyapunov exponent 
usin g th e standard procedure described in detail in 
Ref. [2J|. The tangent space equations are, 

^ + U • V(5Vj + 5U ■ VV' = e6ijj - (V^ + l)^Sip 
ot 

+5N[iP,5ilj], (4) 
V^5C = .g,„ [V (V2,5V') X Vip 

+ V (V^) X V5if\ ■ z, (5) 



s-'('^- 



where 5\5 = V x {5C,i). The nonlinear term for the GSH 
equation is 

5N[^j,5^l?\ = -Z^?5il}, (6) 

and for Manneville's model it is 

5N[iP, ^v] = - (sv-^ + I VV'P) 5V - 2 [(VV') • (v^Vj)] V-- 

(7) 
The magnitude of 5^ is renormalized after a time Im to 
yield a measure of its growth ||5'(/;(tAr)|| which is used to 
calculate the instantaneous Lyapunov exponent, i.e.. 



Ai 



1 

In 



\n\\54>itN)\\. 



(8) 



We have used In — 2 in our simulations. This normal- 
ization is repeated in time to generate many values of 
the instantaneous Lyapunov exponent whose time aver- 
age yields the finite-time Lyapunov exponent 



1 ^' 



Nt 



(9) 



where Nt is the number of renormalizations performed. 
The limit iVf — >■ oo yields the infinite-time Lyapunov ex- 
ponent. 



III. RESULTS 

Experimental measurements of spiral defect chaos have 
typically been performed using large aspect ratio do- 
mains, moderate Rayleigh numbers, and compressed 
gases with Prandtl number of approximately unity [5J- 
Generalizations of the Swift-Hohenberg model, as given 
by Eqs. dl])-©, have been used to study fundamental 
features of spiral defect chaos. The choice of the system 
parameters in the Swift-Hohenberg-type models are im- 
portant in order to yield dynamics that resemble spiral 
defect chaos. In order to estimate the appropriate pa- 
rameters, Xi et al. [9] compared a three-mode amplitude 
equation for the GSH equation with the experimental re- 
sults of Ref. [2^. Using this approach yielded values of 
a = 1, e = 0.7, and gm — 50. It was also chosen to use 
c^ = 2 although a specific physical reason for this par- 
ticular choice is not given. This parameter set has been 



adopted in most numerical explorations of spiral defect 
chaos using the Swift-Hohenberg model, including that 
of Schmitz et al. [Ifl] where it was suggested that spi- 
ral defect chaos in the numerical simulations was only a 
transient. 

The choice of these system parameters highly affects 
the magnitude and dynamics of the mean fiow. The 
magnitude of the mean flow is inversely proportional to 
a and increases with increasing values of the coupling 
strength gm- Cross [26| has presented a careful study of 
the variation of the dynamics with gm- For small val- 
ues of gm the patterns were dominated by target defects 
and for large values of gm the patterns were dominated 
by spirals. Schmitz et al- [lO| computed the appropriate 
value of gm based upon the zig-zag stability boundary 
for e = 0.7 and a — \ and found that gm ~ 12. However 
it was determined by numerical exploration that a larger 
value was required to yield dynamics that resembled spi- 
ral defect chaos and a value of gm = 50 was used in their 
numerics. 

Our calculations here indicate an important role played 
by the parameter c^ on the strength of mean flows. As 
described above in Sec. IIIAl c^ is related to the choice 
of boundary conditions on the bottom and top plates of 
a three-dimensional convection system. Its value physi- 
cally accounts for the viscous damping that occurs near 
the horizontal surfaces. A value of c^ = corresponds to 
perfect slip on the top and bottom plates and a value of 
c^ 7^ corresponds to a no-slip boundary condition. In 
the development of the GSH equation [13| the term c^ is 
introduced as an unknown constant. In the derivation by 
Manneville [IJ, ^M, (? emerges as part of the expansion 
and averaging procedures used when starting from the 
Boussinesq equations. Using the nondimensional form 
of the equations shown in Eq. ([T])-© yields a value of 
(? — 1.03 (after rescaling) for no-slip boundaries [13|- 
However, as pointed out in Ref. [l3|, the precise numeri- 
cal value of c^ depends upon the approximation process 
as well as the manner in which the averaging is done in 
the vertical direction. 

In the following we have explored the details of mean 
flow and spiral defect state using both Manneville's 
model and the GSH equation. What separates our work 
from previous efforts is that we have explored the role of 
the mean flow by systematically varying the value of c^ 
while also computing the leading order Lyapunov expo- 
nent for very long-time simulations. When using Man- 
neville's model we have chosen the system parameters 
to most closely align with those of a Rayleigh-Benard 
convection domain exhibiting spiral defect chaos (i.e., 
e = 0.7, a ~ 2,gm = 50). When using the GSH equation 
we have used the values commonly used in the literature 
(i.e., e = 0.7,(7 = l,gm — 50). For both models we have 
explored the range of c^ values from 0.1 < c^ < 4. We 
note that for our choice of the system parameters the 
Prandtl number is different between our simulations us- 
ing Manneville's model and the GSH model. Our inten- 
tion is not to provide a quantitative comparison between 



these two models but to explore the role of mean flow 
upon spiral defect chaos. 

The variation of the mean flow magnitude with (? can 
be quantified by computing the average kinetic energy 
of the mean flow field. The time dependent value of the 
spatially averaged kinetic energy is given by 



K{t) = \{[Vl^lJl)) 



(10) 



where (t/^:, Vy) are the x and y components of the mean 
flow velocity and (...) represents the spatial average. In 
Fig. [T] we plot {K)^ which represents the time averaged 
value of K(t) as a function of (? . The time averaging is 
performed from the data of the final 10^ time units. The 
error bars represent the standard deviation of the varia- 
tion of (i^)( about its mean value. The circle symbols are 
numerical results for Manneville's model and the square 
symbols are results for the GSH equation. For both mod- 
els the data can be fitted to (K)i oc c~^, indicating the 
rapid increase in the magnitude of the mean flow for de- 
creasing values of c? . 



^ 




0.05 



2 3 

^2 



FIG. 1: The variation of the spatial and time averaged kinetic 
energy iK)^ with (? . The circle symbols are for Manneville's 
model with parameters e = 0.7, (7 = 2, (/,„ = 50, and the 
square symbols are for the GSH equation with e = 0.7, ct = 1, 
g-m ~ 50. The lines are curve fits given by {K)^ = 0.013c~^ for 
the GSH model and {K)^ = 0.005c~^ for Manneville's model. 
The results are averaged using the last 10^ time units of a 
simulation with total time t = 10®. The error bars are the 
standard deviation of the fluctuations of {K)^ about its mean 
value; their maximum values are ~ 10~^ and ~ 10~^ for the 
GSH model and Manneville's model, respectively. 

Figure [2] illustrates the pattern evolution for c^ = 2 
using Manneville's model. The value of c^ = 2 is the 
typical value used in the literature and corresponds to 
a weak mean flow. At small time t — 10'^ the pattern 
is quite complex and contains many dynamic spiral and 
defect structures, similar to the scenario given in the sim- 
ulations of Xi el al. Q . The coarsening to target struc- 
tures is evidenced by the pattern at t = 5 x 10"* and 



represents the approximate duration of the simulations 
by Schmitz et al. using the GSH equation [10]. This 
coarsening process is extremely slow as can be seen by 
comparing the patterns at longer times in Fig. [5] (the 
pattern at i = 10^ for these simulation parameters is 
shown in the bottom right panel of Fig. [3]) . Note that 
during the time evolution some big targets can break up 
and form small defects which will then interact and re- 
combine with other targets or spirals, resulting in the 
process of coarsening. 




t = 1000 



t = 50, 000 




t = 100, 000 



t = 200, 000 



FIG. 2: Instances of the xp field pattern for Manneville's model 
with system parameters e = 0.7, o" = 2, Qm = 50, and c^ — 2, 
showing a slow process of pattern coarsening. Small spiral 
defects dominate at early times (e.g., t — 10"^), whereas at 
large times the pattern has evolved into a state dominated by 
large and very slow moving target structures. 

We have explored the long-time dynamics of the pat- 
terns by varying the strength of the mean flow over the 
parameter range of 0.1 < c^ < 4. Figure [3] illustrates the 
patterns at i = 10^ for four different values of c^. For 
the cases shown, the patterns resemble the state of spi- 
ral defect chaos for c^ < 0.7, showing spatially complex 
structures with rapid dynamics of small-scale spirals and 
localized defects. Although during the evolution process 
some larger spirals or targets may be formed, they are 
transients and will soon break up, with new small spiral 
defects recreated; this procedure will repeat, but no long- 
lasting, coarsened big targets or spirals can exist. How- 
ever, for larger c^ (e.g., — 2) the pattern has coarsened 
to a state dominated by slowly moving target defects, as 
the case given in Fig. [2l These observations indicate 
that strong enough mean flows, as achieved via control- 
ling parameter c^ in the model equations, are needed to 
reach a persistent chaotic or dynamic state, not only for 



breaking up large targets or spirals but also for locally 
recreating new small spiral-type defects to prevent the 
coarsening procedure. 





= 0.1 



c' = 0.5 





c^ = 0.7 d- = 2 

FIG. 3: The long-time patterns for Manneville's model over a 
range of mean flow strengths as determined by the value of (? . 
The simulations are initiated with random initial conditions 
and are integrated until t — 10®. Other parameters are the 
same as Fig. [J] 

To further characterize the properties of system dy- 
namics, in Fig. |4] we examine the time variation of two 
global quantities: the average kinetic energy K(t) as de- 
fined in Eq. (|10|) . and the convective heat flux J(i) which 
is given by 

J(i) = {V) ■ (11) 

Fig. m^a) indicates that as the strength of the mean flow 
increases (i.e., the values of c^ decrease), the magnitude 
of the heat flux decreases with increasing fluctuations 
about the mean value. Such fluctuations are a result of 
the pattern dynamics. More specifically, the dynamics 
of the defects yield the excursions to larger and smaller 
values of the heat flux. For the average kinetic energy 
as shown in Fig. lUJb), both the magnitude of K{f) and 
its fluctuations increase with increasing strength of the 
mean flow. Overall, Fig. U can be used to shed some 
qualitative insight upon the complexity of the dynamics. 
It is clear that for (? = 2 (which is the parameter used in 
most previous studies) the disorder in the kinetic energy 
and heat flux is significantly reduced. 

Although global diagnostics such as the kinetic energy 
and heat flux provide qualitative insights into the com- 
plexity of system dynamics, they can not be used to de- 
termine if the dynamics are chaotic. In order to deter- 
mine if the dynamics are chaotic we have also calculated 
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FIG. 4: Time variation of the pattern dynamics in Man- 
neville's model for different magnitudes of the mean flow. The 
parameters used are the same as those in Fig. |3l (a) Time 
variation of the convective heat flux J(f) given by Eq. (|11[) : 
(b) Time variation of the average kinetic energy Kif) defined 
by Eq. (fTO)l . 



the finite-time leading order Lyapunov exponent using 
the approach discussed in Sec. IIIBI In Fig. [5] we show 
the variation of Ai with the magnitude of the mean flow, 
for which Ai > indicates chaos. The circle symbols are 
the results for Manneville's model and the solid line is 
included only to guide the eye. We emphasize that these 
results are obtained from very long-time simulations to 
ensure that any slow coarsening dynamics have been cap- 
tured. The reported values of Ai are computed using the 
final 10^ time units. The standard deviation of the fluc- 
tuations of Ai about its mean value is included as error 
bars. The maximum value of the error bar for the pa- 
rameters explored is ^ 10^'^. We have also performed 
tests upon the data where Ai is computed over succes- 
sive windows of time of various durations to quantify the 
presence of any slow trends. Our reported values of Ai 
appear to be the converged result and for the duration 
of the simulations explored we did not find any indica- 
tion of slow trends. However, this does not rule out the 
possibility of even slower dynamics not captured in our 
simulations. 



Our simulations indicate that the dynamics are weakly 



chaotic for 



2 where Ai w 0.05. This weak chaos is 



due to the slow dynamics of the pattern, as can be seen 
in Fig. [3] where the targets are slowly moving among the 
more rapid dynamics of defects such as spirals, disloca- 
tions, etc. The value of Ai increases with the decreasing 
value of (? and thus the increasing magnitude of the mean 
flow as given in Fig. [SJ 

We have also included in Fig. [5] the numerical results 
of the GSH equation to enable a more direct comparison 
with prior results in the literature. These results are 
shown as the square symbols and the dashed line. We 
have allowed (? to vary while the remaining parameters 
are the typical values used in the previous studies: e = 
0.7, CT = 1, and gm = 50. Overall our results indicate 
a very similar trend: the dynamics are weakly chaotic 
for c^ w 2 and become increasingly chaotic for smaller 
(? and larger magnitudes of the mean flow. We have 
performed the same series of tests as for Manneville's 
model to account for the presence of any slow coarsening 
process. The maximum magnitude of Ai error bars for 
the GSH equation is ^ 10"^. 



1.2F 




FIG. 5: The variation of the leading order Lyapunov exponent 
Ai with the magnitude of the mean flow. The value of Ai is 
computed using results for the final 10® time units. The circle 
symbols and the solid line are for Manneville's model using 
the system parameters e = 0.7, a = 2, and gm = 50. The 
square symbols and the dashed line are for the GSH equation 
using the typical system parameters of most previous studies: 
e = 0.7, a — 1, and gm — 50. The lines are only meant to 
guide the eye. The error bars represent the standard deviation 
of Ai about its mean value for each value of c^. The maximum 
value of the error bar for Manneville's model is ~ 10~^ and 
for the GSH equation is ~ 10~^. 

Although the connection between the Swift- Hohenberg 
equations and the Boussinesq equations of Rayleigh- 
Benard is phenomenological, a direct comparison be- 
tween the two can provide further physical insights. In 
this comparison we mainly focus on spatial features of 
the mean flow field for chaotic patterns containing a large 
number of spiral defects. In Fig.[BJa) we illustrate the re- 



lationship between the pattern and the magnitude of the 
mean flow. The image shown is for Manneville's model 
at t = 10^ with c^ = 0.1. The dynamics are chaotic with 
Ai w 1.2 as given in Fig. [S] The color contours repre- 
sent the magnitude of the mean flow field |U| where red 
indicates regions of large mean flow and blue indicates 
regions of small mean flow. The roll pattern is shown 
by the black lines which are given by contours of '0 = 0. 
The mean flow tends to reach its local maximum at lo- 
cations that contain defect structures and remains large 
on a length scale of several roll wavelengths around the 
defect. 

We also performed numerical simulations of the three- 
dimensional time-dependent Boussinesq equations that 
describe Rayleigh-Benard convection. We used a paral- 
lel spectral approach that is discussed in detail elsewhere 
(c.f. Q)- We have chosen the reduced Rayleigh number 
e = 0.7, the Prandtl number Pr=l, and an aspect ratio 
of r = 128 for the spatial domain. The top and bottom 
boundaries have the no-slip fluid boundary condition and 
are held at constant temperature, while periodic bound- 
ary conditions are used on the sidewalls. The spatial vari- 
ation of the mean flow is shown in Fig. [Hlb). The mean 
flow is computed as the vertical average of the horizontal 
components of the fluid flow field, with red regions indi- 
cating large mean flow and blue regions corresponding to 
small mean flow. The black lines indicate the contours of 
the convection rolls. Due to the large aspect ratio of the 
domain the computational cost of the numerical simula- 
tion is considerable, and the image shown is for a time 
t — 232.9 where t has been nondimensionalized in the 
usual manner using the time for heat to diffuse across 
the depth of the layer. Qualitatively, the convective roll 
pattern is quite similar to what is shown for the Swift- 
Hohenberg-type equations. The spatial variation of the 
mean flow fleld, on the other hand, is different in terms of 
its rate of decay with the distance from a defect core. The 
mean flow fleld is largest at the core regions of the defect 
structures and decays rapidly with distance away from 
the defect core. In most cases the length scale of this de- 
cay is approximately that of a single roll (see Fig. Eljb)), 
as compared to several rolls for results of Manneville's 
model shown in Fig. [Hlja). 

The derivation of the Swift-Hohenbcrg equations is 
based upon a long wavelength approximation and hence 
the corresponding results are not expected to be accurate 
in the core regions of defect structures Q • This is illus- 
trated in Fig. [7] which shows a close-up view of the mean 
flow structure near a spiral defect, with panel (a) for re- 
sults of Manneville's model and panel (b) for the Boussi- 
nesq equations. In both cases the spiral is rotating in 
a counter-clockwise direction. For the Swift-Hohenberg 
equations the mean flow exhibits a quadrupole struc- 
ture centered upon the spiral. The magnitude of the 
quadrupole is spatially asymmetric and varies with the 
defect dynamics. The mean flow for the Boussinesq equa- 
tions is a vortex rotating in the opposite direction to the 
rotation of the spiral. Our results for the Boussinesq 




FIG. 6: (color) Spatial variation of the mean flow field for 
chaotic patterns containing many spiral defects. The color 
contours represent the magnitude of the mean flow field, with 
red regions corresponding to large mean fiow magnitude and 
blue regions corresponding to small mean flow magnitude. 
The black contours indicate convective roll boundaries, (a) 
Results for Manneville's equations at i = 10® using system 
parameters (e = 0.7, a — 2, gm = 50, c^ — 0.1). (b) Results 
for the Boussinesq equations using system parameters e = 0.7, 
Pr= 1, and F = 128 at time t = 232.9. 



equations are in agreement with those presented by Bo- 
denschatz et. al (see Fig. 18 in Ref. Q). 



IV. CONCLUSION 

Using generalizations of the Swift-Hohenberg equation 
we have explored the spiral defect chaos state and the role 




FIG. 7: (color) A close-up view of the spatial variation of 
the mean flow near a spiral defect in a chaotic pattern. The 
color contours represent the magnitude of the mean flow as in 
Fig. [6l the black contours indicate the convective roll bound- 
aries, and the arrows are vectors for the mean flow, (a) Re- 
sults using Manneville's equations with system parameters 
e = 0.7,0- = 2, g„ = 50, and c^ = 0.4 at time f = 7 x 10^. 
(b) Results from the Boussinesq equations using system pa- 
rameters e = 0.7, a — 1, and F = 128 at time t = 232.9. 



of the mean flow. Our results show that it is possible to 
generate chaotic dynamics using Swift-Hohenberg-type 
model equations that resemble the spiral defect chaos of 
Rayleigh-Benard convection. The important insight is 
that the strength of the mean flow must be large enough. 
Reasonable parameter values of e, a, and gm that lead 
to spiral defect chaos have been explored at some length 
in previous literature. However, the precise value of pa- 
rameter c^ depends on the approximations used in the 
derivation of Manneville's equations and the appearance 
of c^ is phenomenological in writing down the GSH equa- 



tions. Our results show that the dynamics vary strongly 
with the magnitude of (? and the commonly used value 
oi (? — 2 yields a mean flow that is not strong enough to 
generate persistent dynamics that resemble spiral defect 
chaos. By choosing a smaller value of (? the dynamics 
are chaotic for as long as we have simulated (t — 10^ time 
units). Although we focused our discussion upon results 
generated using Manneville's equations, our conclusions 
and insights also apply to the GSH equations. The par- 
ticular choice of the form of the nonlinearity does not 
strongly affect the results in any significant way that we 
have found. The necessity of a strong enough mean flow 
to support spiral defect chaos is in agreement with what 
has been found from the Boussinesq equations [l3] . 

However, there are significant differences between nu- 
merical results from the Swift-Hohenberg-type model 
equations and those generated using the full three- 
dimensional Boussinesq equations. The Swift-Hohenberg 
equations are not expected to capture the small scale fea- 
tures correctly due to the long- wavelength approximation 
involved, and we have quantified some aspects of this by 
comparing the mean flow fields around a single spiral de- 



fect. 

We anticipate that our results will have several uses. 
First is that Swift-Hohenberg-type equations can be used 
as a model to study spatiotemporal chaos in a sys- 
tem with direct relevance to fiuid phenomena such as 
Rayleigh-Benard convection. This is a significant ad- 
vantage since numerical simulations of the Boussinesq 
equations are computationally very expensive J\. The 
Swift-Hohenberg model equations could be used to ex- 
plore fundamental ideas of spatiotemporal chaos that are 
currently inaccessible to the full fluid equations. Exam- 
ples include a detailed study of microextensivity [23, |2M| 
or the cornputation of the characteristic Lyapunov expo- 
nents [29|, |30| . In addition, our comparison of the spatial 
variation of the mean flow around spiral defect core struc- 
tures could be used to guide the development of more 
accurate theoretical descriptions of spiral defect chaos. 
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